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Abstract 

Coffee is one of the world's most traded agricultural products. Modeling studies have predicted that climate change 
will have a strong impact on the suitability of current cultivation areas, but these studies have not anticipated possible 
mitigating effects of the elevated atmospheric [C0 2 ] because no information exists for the coffee plant. Potted plants 
from two genotypes of Coffea arabica and one of C. canephora were grown under controlled conditions of irradiance 
(800 umol nr 2 s -1 ), RH (75%) and 380 or 700 uL C0 2 L 1 for 1 year, without water, nutrient or root development 
restrictions. In all genotypes, the high [C0 2 ] treatment promoted opposite trends for stomatal density and size, which 
decreased and increased, respectively. Regardless of the genotype or the growth [C0 2 ], the net rate of C0 2 
assimilation increased (34-49%) when measured at 700 than at 380 uL C0 2 L 1 . This result, together with the almost 
unchanged stomatal conductance, led to an instantaneous water use efficiency increase. The results also showed a 
reinforcement of photosynthetic (and respiratory) components, namely thylakoid electron transport and the activities 
of RuBisCo, ribulose 5-phosphate kinase, malate dehydrogenase and pyruvate kinase, what may have contributed to 
the enhancements in the maximum rates of electron transport, carboxylation and photosynthetic capacity under 
elevated [COJ, although these responses were genotype dependent. The photosystem II efficiency, energy driven to 
photochemical events, non-structural carbohydrates, photosynthetic pigment and membrane permeability did not 
respond to [COJ supply. Some alterations in total fatty acid content and the unsaturation level of the chloroplast 
membranes were noted but, apparently, did not affect photosynthetic functioning. Despite some differences among 
the genotypes, no clear species-dependent responses to elevated [COJ were observed. Overall, as no apparent sign 
of photosynthetic down-regulation was found, our data suggest that Coffea spp. plants may successfully cope with 
high [COJ under the present experimental conditions. 
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Introduction 

Over the past 800,000 years, atmospheric [COJ has varied 
between 180 uL C0 2 L" 1 (glacial periods) and 280 uL C0 2 L 1 
(interglacial periods) as Earth moved between ice ages. From 
pre-industrial levels of 280 pL C0 2 L" 1 , [COJ has increased 
steadily to 384 uL C0 2 L" 1 in 2009 (reached 400 uL C0 2 L \ 
measured in Mauna Loa Observatory in Hawaii in 2013), and 



levels are estimated to reach between 450 and 600 pL L" 1 by 
the year 2050 and between 730 and 1020 pL L 1 by 2100, 
depending on future scenarios of anthropogenic emissions 
[1.2]. 

Changes in atmospheric [COJ affect fundamental plant 
processes and may alter plant growth, agronomic yields and 
quality [3-6]. Crops sense and respond directly to rising 
atmospheric [COJ through changes in photosynthesis and 
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stomatal conductance (g s ), which is the basis of what can be 
considered a C0 2 fertilization effect on crop yield [7]. Such an 
effect may eventually strengthen the plant and alter its 
tolerance limits to environmental constraints. For example, it is 
known that high temperatures can reduce net C gain by 
increasing photorespiration. Additionally, a high [C0 2 ] supply 
can increase the net photosynthetic rate (P n ) in C 3 plants (often 
above 50%) due to the higher carboxylation rate of 
ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCo), 
which results from the simultaneous increases in substrate 
availability and competitive inhibition to 0 2 [3,7-10]. Therefore, 
by reducing photorespiration, C0 2 enrichment is expected to 
enhance P n to a greater degree at high than at low 
temperatures, thereby, at least partially, offsetting the effects of 
supra-optimal temperatures on yield [2,1 1,12]. 

Overall, g s is consistently but not universally decreased at 
elevated [C0 2 ] [13-15] what may decrease the transpiration 
rates (T r ). Because increasing atmospheric [C0 2 ] enlarges the 
gradient that ensures adequate diffusion of C0 2 from the 
atmosphere to the chloroplasts, a rise in [C0 2 ] should produce 
greater P n coupled with lower T r , which would ultimately 
improve the water-use efficiency (WUE) in a large number of 
plant species [8,14,16]. Furthermore, decreases in stomatal 
opening, stomatal density (SD) and stomatal index (SI), all of 
which contributing to a reduction in g s , have been reported 
when plants are grown in elevated [C0 2 ] [8,17-19], although in 
some cases g s and not SD determines the long-term reduction 
in leaf T r [20]. Nonetheless, the anticipated stomatal closure 
that is usually observed at elevated [C0 2 ] will inevitably be 
associated with lower latent heat loss, thereby increasing leaf 
temperatures [2]. 

The degree of downregulation (acclimation) of 
photosynthesis in response to C0 2 enrichment is variable 
among plants, depends on the interactions with other 
environmental limitations [7,10,21] and even changes with the 
developmental stage of the same plant [22]. If present, the 
downregulation of photosynthesis decreases (but does not 
completely eliminate) the stimulation effect of elevated [C0 2 ] 
[16]. Concomitantly, C0 2 enrichment may lead to increases in 
N- and water-use efficiency and decreases in leaf and plant 
nutrient concentrations. Yet a strong P n increase under high 
growth [C0 2 ] has often been reported to enhance relative 
growth rate to a much lower extent (approximately 10%). This 
discrepancy may result from a limitation on sink strength that 
prevents the plant from fully utilizing the higher photosynthate 
production due to, e.g., a limitation of meristematic tissue 
related to a deterministic growth pattern. Such a limitation 
could lead to an increase in leaf nonstructural carbohydrates 
(NSC) associated with a lower export rate to other tissues 
[10,23], and implicating a reduction of P : regeneration in the 
chloroplast [5,24]. Such NSC rise could trigger a signaling 
mechanism, promoting a cascade of molecular and 
biochemical responses, e.g., the depression of gene 
expression and the amount/activity of photosynthetic enzymes, 
including that of RuBisCo, or a reduction in the levels of all 
components of the photosynthetic apparatus [21,23]. That 
would contribute to the lower the rates of net assimilation and 
photosynthetic capacity, linked to reductions in the maximum 



apparent carboxylation velocity, V cmax , and the in vivo maximum 
apparent rate of electron transport, J max [5,7-9,19,22,25]. In 
fact, the reduction of photosynthesis has been mechanistically 
and quantitatively attributed to decreased V cmax and investment 
in RuBisCo [13], but the responsible factors may also include 
reduced ribulose 1 ,5-bisphosphate (RuBP) regeneration, which 
decreases J max due to lowered electron transport capacity or P; 
availability in the chloroplast for ATP synthesis [3,8,10,21,22]. 
Trees and shrubs, particularly fast-growing species, generally 
have a larger sink capacity (root-trunk system) than annuals 
and usually show higher photosynthetic stimulation when 
grown at high [C0 2 ] than do shrubs and annual crops [8,13,26]. 

Coffee is a tropical crop that is currently grown in 
approximately 80 countries, making it one of the world's most 
traded agricultural products. The sale of coffee generates over 
US $90,000 million each year and is the economic basis of 
many tropical developing countries. The genus Coffea 
comprises more than 100 species, among which C. arabica L. 
and C. canephora Pierre ex A. Froehner that, together, are 
responsible for approximately 99% of world coffee bean 
production [27,28]. Due to the world's ongoing climatic 
changes, there has been increasing concern regarding the 
suitability of traditional coffee producing areas. Modeling 
studies of predicted climate scenarios, mostly related to 
increased air temperatures, have estimated dramatic effects on 
this crop, including severe yield losses in Mexico [29], 
extensive reductions of suitable areas in Brazil [30] and the 
extinction of wild populations of C. arabica in Ethiopia [31]. 
However, these negative predictions of the effects of global 
climate changes have not considered the mitigating effects of 
increasing atmospheric [COJ on the harmful impacts of 
elevated temperatures on the crop. Moreover, the ability that 
coffee plants often display to successfully adjust their 
metabolism to stressful environmental conditions [27,32-34] 
must also be taken into account. 

Despite the agronomic importance of coffee, no information 
exists regarding the effects of C0 2 on the physiology of this 
crop. The coffee plant displays low P n , typically in the range of 
4-11 umol m 2 s 1 , at current atmospheric [COJ and saturating 
light conditions [2], which is largely due to diffusive, rather than 
biochemical, limitations to photosynthesis [35,36]. It has also 
been proposed that coffee plants are, within given limits, 
capable of avoiding the downregulation of photosynthesis 
through their high capacity for starch accumulation [36,37], and 
thus we hypothesize that coffee will sustain relatively high P n in 
a scenario of increasing atmospheric [COJ. Here, we report the 
first results concerning the underlying mechanisms associated 
with the responses of the photosynthetic apparatus to elevated 
atmospheric [COJ in coffee. We investigated three important 
genotypes from C. arabica and C. canephora grown at either 
380 or 700 pL C0 2 L 1 air under controlled conditions without 
water, nutrient or root development restrictions. Our results are 
discussed in the context of current models of photosynthetic 
performance in a scenario of increasing atmospheric [COJ. 



PLOS ONE | www.plosone.org 



2 



December 2013 | Volume 8 | Issue 12 | e82712 



Enhanced [C0 2 ] Impact on Coffea spp. 



Materials and Methods 

Plant material and experimental design 

Potted plants of 1 .5 years in age from C. arabica L. (cv. Icatu 
and cv. Catucai IPR 108) and C. canephora Pierre ex Froehner 
cv. Conilon Clone 153 (CL 153), grown in 12 L pots, were 
transferred from a greenhouse (ambient [C0 2 ]) into walk-in 
growth chambers (EHHF 10000, ARALAB, Portugal). The 
plants were then grown under controlled environmental 
conditions of temperature (25/20°C, day/night), RH (75%), 
irradiance (approximately 800 pmol m 2 s 1 ), photoperiod (11.5 
hours) and two C0 2 concentrations (380 or 700 uL L 1 ). Each of 
the 10 plants per genotype and [C0 2 ] treatment was fed on a 
monthly basis with 5 g of the following fertilizer mixture: 7% 
Ca(N0 3 ) 2 , 5% KN0 3 , 7.8% P 2 0 5 , 17% K 2 0, 1.6% MgO, 20% 
MgS0 4 , 0,02% H3BO3 and 0.01% ZnS0 4 . To reinforce the N 
and Ca availability, a complementary fertilization of 2 g was 
conducted every 3 months with a mixture of 27% NH 4 N0 3 and 
6% CaO. Both fertilizers were provided as solid spheres that 
slowly dissolved over successive watering, allowing a gradual 
release of minerals to the soil/plant. To complement the 
availability of micronutrients, 500 mL of a solution containing 
0.02% Fe-EDTA, 0.01% CuS0 4 , 0.01% MnCI 2 , and 0.005% 
H 2 Mo0 2 , were added on a monthly basis. 

After one year, a large set of parameters was evaluated in 
recently mature leaves from the top branches (light exposed) 
that had fully developed under each [COJ. For biochemical 
analysis, leaf material was collected after approximately 2 
hours of illumination from 6 to 8 plants of each genotype, flash 
frozen in liquid nitrogen and stored at -80°C until analysis. 
Whenever possible, all analyses were performed on the same 
leaves. There was no apparent restriction on root development, 
as judged by visual examination at the end of the experiment 
when plants were removed from their pots. 

Stomatal determinations and specific leaf area 

As stomatal traits (e.g., stomatal density) are responsive to a 
range of environmental factors in addition to their relative 
position on the leaf [38], observations were made on similar 
areas of the leaf lamina by taking imprints from the abaxial leaf 
surface at the point of maximal leaf width, from the margin to 
the main central leaf vein, using colorless nail polish and 
adhesive transparent cellophane tape. The imprints were then 
placed on glass slides and examined under a light microscope 
(Olympus BX50, Japan) at 400x magnification, attached to a 
digital camera (Camera Zeiss Axiovision, Germany). Five 
samples per replicate per genotype, with three fields of view 
within each sampled area, were analyzed. 

The stomatal traits were analyzed following [39]. Stomatal 
density (SD) was calculated as the number of stomata per leaf 
area unit, whereas the stomatal index (SI) was calculated as SI 
= [(stomata )/(total cells +stomata)] *100. For stomatal size 
(SS), 30 randomly selected stomata were measured using an 
ocular micrometer, and their areas were calculated as SS = 
Trab, where a and b are 1/2 length and 1/2 width, respectively, 
thus assuming that the stomatal shape is an ellipse. 

The specific leaf area (SLA) was determined in 5 samples of 
10 leaf discs (0.5 cm 2 each) after drying at 80°C for 24 h. 



Gas exchanges 

The responses of the net photosynthetic rate (P n ) to internal 
C0 2 concentration (C;) (P n /Ci curves) were assessed using a 
portable open-system infrared gas analyzer (Li-Cor 6400, 
LiCor, Lincoln, USA). The curves were obtained from 5 to 7 
plants per genotype under the light conditions of the growth 
chambers (ca. 800 pmol m 2 s 1 ) using a 10-12 stepwise 
external [COJ levels from 50 to 1800 pL L 1 . Each step value 
was taken only when P n , stomatal conductance to water vapor 
(g s ) and transpiration rate (T r ) were stable, which was 
approximately 5-7 min after each [COJ level was imposed. 
From the P n /Ci curves, were estimated the in vivo maximum 
apparent rate of carboxylation (V cmax ), the in vivo maximum 
apparent rate of electron transport (J max ), the respiration in the 
presence of light (R d ) and the triose phosphate utilization rate 
(TPU), as described elsewhere [40,41]. As the conditions 
inside the growth chamber were quite stable (including air 
humidity kept at 75%) and no limitations (e.g., water, nutrients) 
were imposed, the leaf temperature was no higher than 1 °C 
above air temperature (25 °C). Therefore, in the experimental 
conditions, the leaf-to-air vapor pressure deficit was also quite 
constant and low (the highest value was close to 0.985 kPa). 

The steady-state values of leaf P n , g s and T r were taken from 
the same P n /C; curves at 380 or 700 pL C0 2 L 1 . Leaf 
instantaneous water-use efficiency (iWUE) was calculated as 
the P n -to-T r ratio, representing the units of assimilated C0 2 per 
unit of water lost through transpiration. 

The photosynthetic capacity, A max (representing the light- and 
C0 2 -saturated rate of photosynthesis under optimal 
temperature), was measured through 0 2 evolution in a Clark- 
type 0 2 electrode (LD2/2, Hansatech, UK) using leaf discs 
(1 .86 cm 2 ), following the procedures described elsewhere [42]. 
A max was obtained at 25 °C under saturating [COJ conditions 
(ca. 7%, supplied by 400 pi KHC0 3 , 2 M) by exposing the leaf 
samples to increasing irradiances up to 1500 pmol m 2 s 1 using 
a Bjorkman lamp (Hansatech) and neutral filters. 

Chlorophyll a fluorescence analysis 

Chlorophyll (Chi) a fluorescence parameters were 
determined on the same leaves used for the gas exchange 
measurements using a PAM-2000 system (H. Walz, Effeltrich, 
Germany), as previously described [42], following formulae 
discussed elsewhere [43,44]. Briefly, the estimation of the 
maximal photochemical efficiency of photosystem (PS) II 
(F v /F m ) was performed on overnight dark-adapted leaves, using 
a 0.8 s saturating light pulse of 7500 pmol m 2 s _1 . The 
photochemical quenching (q P ), which denotes the proportion of 
energy trapped by PSII and driven to photochemical events, 
and the maximal PSII efficiency of energy conversion under 
light (F v 7F m ") were determined under photosynthetic steady- 
state conditions, at the light conditions of the growth chambers, 
using superimposed saturating flashes. 

Thylakoid electron transport rates 

Subchloroplast fractions were obtained from a pool of leaf 
material (ca. 5 g FW) from 5 or 6 plants and processed as 
previously optimized for coffee leaves [33,45]. The in vivo 
electron transport rates associated with both PSI 
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(DCPIPH 2 ^MV) and PSII, including (H 2 O^DCPIP) or 
excluding (DPC— >DCPIP) the oxygen evolving complex (OEC), 
were measured polarographically using an LW2 0 2 electrode 
(Hansatech) at 25 °C. The assays were performed using 1 mL 
of the reaction mixture (containing ca. 100 mg Chi) and a 
PPFD of approximately 3000 pmol m 2 s 1 given by a Bjorkman 
lamp (Hansatech). 

Enzyme activities 

Four freshly cut leaf discs (0.5 cm 2 each) were used to 
measure the activity of several enzymes involved in carbon 
metabolism. Each sample was homogenized in a cooled mortar 
using 100 mg insoluble PVPP and 1 mL of the extraction buffer 
100 mM Tris-HCI (pH 8), which contained 10 mM MgCI 2 , 10 
mM NaHC0 3 , 10 mM fi-mercaptoethanol, 2 mM DTT, 1% (v/v) 
Triton X-100, 10% (v/v) glycerol and a "complete-protease 
inhibitor cocktail" 2% (v/v) designed to protect the enzymes 
from protease action (Roche, ref. 04693159001). The extracts 
were centrifuged (16,000 g, 20 min, 4 °C) and the supernatant 
was used for the enzyme assays, all of which were based on 
NADH oxidation at 340 nm, at 25 °C, in 1 mL final volume in 
the cuvette. 

The total ribulose-1 ,5-bisphosphate carboxylase/oxygenase 
(RuBisCo: EC 4.1.1.39) activity was determined following 
previously described procedures [46], with some modifications: 
20 pL of the leaf extract were added to a cuvette containing an 
assay medium of 50 mM Tris-HCI, pH 8.0, 15 mM MgCI 2 , 20 
mM NaHC0 3 , 100 mM creatine phosphate, 10 mM ATP, 0.2 
mM NADH, 20 U mL 1 creatine kinase, 15 U mL 1 3- 
phosphoglycerate kinase and 15 U mL 1 glyceraldehyde-3- 
phosphate dehydrogenase. After an incubation period of 10 
min, 50 uL of 20 mM RuBP were added to initiate the reaction. 

The activity of ribulose 5-phosphate kinase (RuB5PK: EC 
2.7.1.19) was determined following procedures described 
elsewhere [41], by adding 20 pL of the leaf extract to a cuvette 
containing an assay medium of 100 mM Tris-HCI, pH 8.0, 8 
mM MgCI 2 , 40 mM KCI, 20 mM phosphoenolpyruvate, 5 mM 
ATP, 1 mM NADH, 20 mM DTT, 8 units pyruvate kinase, 10 U 
mL 1 lactate dehydrogenase and 5 U mL 1 
phosphoriboisomerase. After an incubation period of 15 min, 
the reaction was initiated by adding 10 pL of 500 mM ribose-5- 
phosphate. 

The pyruvate kinase (PK: EC 2.7.1.40) activity was 
determined as previously described [47], with some changes: 
20 pL of the leaf extract were added to a cuvette containing an 
assay medium of 100 mM Tris-HCI (pH 7.0), 10 mM MgCI 2 , 0.2 
mM NADH, 1 mM fructose-1,6-bisphosphate, 45 mM ADP and 
6.3 U mL 1 lactate dehydrogenase. The reaction was initiated 
by adding 100 pL of 10 mM phosphoenolpyruvate. 

The determination of NADH-dependent malate 
dehydrogenase (MDH: EC 1.1.1.37) activity followed 
procedures described elsewhere [48], with some changes, by 
adding 20 pL of the leaf extract to a cuvette containing an 
assay medium of 50 mM Tris-HCI (pH 8.0) and 0.1 mM NADH. 
The reaction was initiated by adding 20 pL of 20 mM 
oxaloacetate as substrate. 



Non-structural carbohydrate quantification 

Soluble sugars were determined in approximately 150 mg of 
powdered frozen material, based on the method previously 
described [49]. The samples were homogenized in 2 mL of cold 
H 2 0, left to extract for 20 min on ice and centrifuged (12,000 g, 
5 min, 4 °C). The supernatant was boiled to denature the 
proteins (3 min), placed on ice (6 min) and centrifuged again. 
The obtained clear solution was then filtered (0.45 urn, nylon) 
before the injection of a 50 pL aliquot into an HPLC system 
equipped with a refractive index detector (Model 2414, Waters, 
USA). The separation of sugars was performed using a Sugar- 
Pak 1 column (300 x 6.5 mm, Waters) at 90 °C, with H 2 0 as 
the eluent (containing 50 mg EDTA-Ca L 1 H 2 0) and a flow rate 
of 0.5 mL min -1 . To overcome the presence of non-pure peaks 
from this separation, another 20 pL aliquot of each sample was 
injected through a DionexCarboPac PA1 analytical column (4 x 
250 mm, Thermo Scientific, USA) coupled to a 
DionexCarboPac PA1 Guard (4 * 50 mm) at 20 X. Ultrapure 
water and 300 mM NaOH were used as eluents (water from 0 
to 50 min; NaOH from 50 to 65 min; and water from 65 to 80 
min for re-equilibration) at a 1 mL min 1 flow rate. Standard 
curves were used for the quantification of each sugar. 

Starch determination was conducted according to previously 
described procedures [50], with some changes. After adding 1 
mL of boiling millipore water to 100 mg of frozen leaf material, 
the samples were immediately placed into boiling water, 
shaken for 10 min and centrifuged (10,000 g, 2 min, 4°C). The 
supernatant containing the soluble sugars was discarded, and 
the insoluble pellet containing starch was collected, further 
washed and centrifuged twice in 1 mL millipore water. The 
insoluble residue was thoroughly homogenized in 1 mL H 2 0, 
and the resulting suspension was autoclaved (120°C, 3 h) to 
promote starch gelatinization. Starch hydrolysis was performed 
in 0.5 mL of the solution by adding 0.5 mL of 400 mM 
citrate/KOH buffer, pH 4.6, containing 30 U amyloglucosidase 
(EC 3.2.1.3, Sigma) and 2 U a-amylase (EC 3.2.1.1, Sigma), 
followed by an overnight shaking incubation at 37 °C. The 
extract was then cleared by centrifugation (10,000 g, 2 min, 
4°C), and the glucose derived from starch was enzymatically 
determined on a 20 pL aliquot of the supernatant, using 100 
mM imidazol buffer, pH 6.9, containing 1.1 mM ATP, 0.5 mM 
NADP, 1 U hexokinase (EC 2.7.1.1) and 1 U glucose-6- 
phosphate dehydrogenase (EC 1.1.1.49) in a 1 mL final 
volume. The determination was conducted 
spectrophotometrically at 340 nm. 

Photosynthetic pigments 

Total chlorophylls and carotenoids were extracted from four 
freshly cut leaf discs (0.5 cm 2 each) using 80% (v/v) aqueous 
acetone and quantified according to procedures described 
elsewhere [51]. 

Membrane permeability 

Ten freshly cut leaf discs (0.5 cm 2 each) were rinsed 3 times 
(1 min) with demineralized water and subsequently floated on 
10 mL of demineralized water at 20°C, following previously 
described procedures [52]. The electrolyte leakage was 
measured until leakage stabilization occurred, at 22 h, using a 
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conductivity meter (Crison GLP31, Crison Instruments, S.A., 
Spain). Total conductivity was obtained after the flasks were 
exposed to 90°C for 2 h in an oven and cooled. Membrane 
leakage was given as a percentage of total conductivity. 

Lipid quantification from chloroplast membranes 

The lipid fraction from the enriched chloroplast membranes 
was obtained from 3 to 4 g (FW) of leaf tissue, followed by the 
quantification and identification of fatty acids (FAs), as has 
been previously described for coffee leaves [28]. The value for 
the total fatty acid (TFA) content corresponds to the sum of 
individual FAs, while the double bond index (DBI) was 
calculated as DBI = [(% monoenes + 2 x % dienes + 3 x % 
trienes / (% saturated FAs)]. 

Statistical analysis 

The various measured and calculated parameters were 
analyzed using two-way ANOVAs (P < 0.05) to evaluate the 
differences between C0 2 treatments or genotypes, followed by 
a Tukey test for mean comparisons among genotypes within 
the same [C0 2 ] condition and an F test for mean comparisons 
between [C0 2 ] treatments within the same genotype; a 95% 
confidence level was adopted for all tests. For the sake of 
simplicity concerning the ANOVA results, in the figure and table 
captions for each parameter are only indicated when significant 
the differences related to growth [C0 2 ], genotype of for the 
[C0 2 ] x genotype interaction. 

Results 

Stomatal traits and specific leaf area 

Preliminary observations confirmed that coffee plants have 
stomata only on the abaxial leaf surface. A consistent trend for 
lower stomatal density, SD (5-14%), and higher stomatal size, 
SS (3-7%), was found under high [C0 2 ] in all genotypes, 
although statistical significance was reached only in the 
genotype Icatu (Table 1 ). Regardless of genotype, the stomatal 
index (SI) and specific leaf area (SLA) did not respond 
significantly to the [C0 2 ] treatments. 

Leaf gas exchanges 

Irrespective of genotype, a tendency to higher values was 
observed for the net photosynthetic rate (P n ) under higher 
[C0 2 ] (between 34% in IPR 108 and 49% in CL 153). 
Furthermore, P n measured at 380 or 700 uL C0 2 L -1 showed 
similar values in the plants grown either at normal or elevated 
[C0 2 ], thus the changes in P n were independent of growth 
[C0 2 ] conditions (data not shown). 

Concerning the stomatal conductance to water vapor (g s ), no 
significant changes were observed between the [C0 2 ] 
treatments for each genotype. However, g s values tended to be 
lower when measured at 700 uL C0 2 L \ from 4% (CL 153) to 
28% (Icatu) (Figure 1). Again, the non-significant decreasing 
trend observed in g s when measured at 700 uL C0 2 L" 1 was 
irrespective of growth [C0 2 ] (data not shown). 

Despite the absence of significant changes in the above 
parameters when measured at 700 uL C0 2 L" 1 in comparison to 



Table 1. The stomatal density (SD), stomatal size (SS), 
stomatal index (SI) and specific leaf area (SLA) in the 
leaves of C. arabica (Icatu and IPR 108) and C. canephora 
(Conilon CL 153) grown under 380 or 700 uL C0 2 L 1 . 



Genotype 


CL 153 




Icatu 




IPR 108 




C0 2 


380 |JL 


700 |JL 


380 uL 


700 uL 


380 uL 


700 uL 


Treatment 


L- 1 


L- 1 


L- 1 


L- 1 


L- 1 


L- 1 


Stomatal 
Density 


195.2 ar 


174.8 ar 


193.9 ar 


165.9 as 


190.7 ar 


180.6 ar 


n. stomata 
mm" 2 ) 


8.8 


6.8 


8.4 


7.2 


4.7 


3.6 


Stomatal 
ndex 


14.0 br 


14.8 br 


24.0 a ' 


20.7 a ' 


26.2 a ' 


22.6 a ' 


(%> 


0.4 


0.4 


1.1 


0.8 


1.3 


0.5 


Stomatal 
Size 


273.2 c ' 


289.8 cr 


349.7 b 5 


374.8 b ' 


398.1 a ' 


409.3 a ' 


(Mm 2 ) 


7.7 


7.5 


7.5 


7.3 


6.1 


9.1 


SLA 


13.5 ar 


13.6 ar 


13.4 ar 


11.6 a r 


15.8 ar 


12.9 ar 


(m 2 Kg" 1 ) 


2.2 


3.8 


1.4 


1.2 


1.8 


1.4 



For each parameter, the mean values ± SE (n = 30 for SD, SS and SI; n = 5 for 
SLA) followed by different letters express significant differences between cultivars 
for the same CO2 treatment {a, b, c) or between CO2 treatments within the same 
cultivar (r, s). The ANOVA for SD showed significant differences only between CO2 
treatments within the same cultivar; that for SI showed significant differences 
between cultivars for the same CO2 treatment; and that for SS showed significant 
differences between cultivars for the same CO2 treatment and between CO2 
treatments within the same cultivar; that for SLA did not show significant 
differences. 

doi: 10.1371/journal.pone.0082712.t001 

380 uL C0 2 L" 1 , the opposite tendencies of g s (and T r ) and P n 
when measured at high [C0 2 ] led to significant increases in the 
instantaneous water use efficiency (iWUE) in all genotypes, 
between 56% and 112% (Figure 1). 

Among the genotypes, differences in the functioning of the 
photosynthetic apparatus were found. Both the in vivo apparent 
rates for the maximum carboxylation (V cmax ) and maximum 
electron transport (J max ) were unresponsive to C0 2 enrichment 
in CL 153 and IPR 108, whereas in Icatu these parameters 
increased 52% and 37%, respectively, at elevated relative to 
normal growth [C0 2 ] (Table 2). The balance between RuBisCo 
carboxylation and electron transport (V cmax /J max ratio) was 
unaltered by growth [C0 2 ], averaging 1.21 across all 
genotypes. Similarly, the triose phosphate utilization for 
sucrose and starch synthesis (TPU), dark respiration rate, R d 
(Table 2) and photosynthetic capacity, A max (Figure 2), were 
also unresponsive to the applied C0 2 treatments, although A max 
tended to increase with C0 2 in CL 153 (18%) and Icatu (25%). 

Chlorophyll a fluorescence analysis 

Only modest changes were noted for the maximal 
photochemical efficiency of PSII (F v /F m ) in the plants grown 
under high [C0 2 ], with a maximal reduction of 3% in Icatu 
(Table 3). These values were not far from those expected 
under shade conditions (0.78), where PSI makes a larger 
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Figure 1. Gas exchanges under ambient [C0 2 ]. Changes in leaf net photosynthesis (P n ), stomatal conductance to water vapor 
(g s ) and instantaneous water use efficiency (iWUE) in C. arabica (Icatu and IPR 108) and C. canephora (Conilon CL 153) measured 
at their growth [C0 2 ]: 380 (white bars) and 700 (black bars) uL C0 2 L 1 . For each parameter, the mean values + SE (n = 6-8) 
followed by different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 
treatments within the same cultivar (r, s). The ANOVA for P n , showed significant differences between cultivars for the same growth 
C0 2 treatment and between C0 2 treatments within the same cultivar; that for iWUE showed significant differences between C0 2 
treatments within the same cultivar; that for g s did not show any significant difference,. 

doi: 10.1371/journal.pone.0082712.g001 



contribution to the F 0 level [53]. Regardless of C0 2 treatment, 
no significant changes were found under photosynthetic 



steady-state conditions for either photochemical quenching (q P ) 
or the photochemical efficiency of PS II (F v 7F m '). 
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Table 2. Variation of the estimations of the maximum rate 
of carboxylation (V cmax ), the maximum rate of carboxylation 
limited by electron transport (J max ), the rate of respiration in 
the presence of light (R d ), and the triose-phosphate (TPU) 
rate of utilization (all calculated from the P n /Ci curves and 
expressed in umol C0 2 nr 2 s -1 ), as well as the J max /V cmax 
ratio, in the leaves of C. arabica (Icatu and IPR 108) and C. 
canephora (Conilon CL 153) grown under 380 and 700 uL 
L-1 of C0 2 . 



Genotype 


CL 153 




Icatu 




IPR 108 




CO2 Treatment 


380 uL 
L" 1 


700 ul_ 
L- 1 


380 uL 
L" 1 


700 uL 

L' 1 


380 uL 
L" 1 


700 uL 
L" 1 


Vcmax 


50.2 a ' 


48.5 ar 


34.8 b 5 


53.0 ar 


43.7 a ' 


51.8 ar 


(|jmol CO2 m" 2 


4.5 


2.1 


2.2 


8.6 


2.1 


5.8 


Jmax 


58.6 a ' 


62.0 a r 


45.1 b 5 


62.0 a ' 


52.7 ab ' 


53.4 a ' 


(Mmol CO2 m" 2 


2.1 


8.5 


2.6 


6.7 


2.5 


6.0 


J max W c max 


1.25 ar 


1.25 ar 


1.31 ar 


1.25 ar 


1.22 ar 


1.12 ar 




0.15 


0.09 


0.10 


0.10 


0.06 


0.10 


Rd 


2.11 ar 


2.38 a ' 


1.74 ar 


3.68 a ' 


1.80 ar 


3.14 ar 


(umol CO2 m" 2 


0.76 


0.45 


0.46 


0.54 


0.41 


0.61 


TPU 


3.19 ar 


3.51 a r 


2.91 a ' 


3.68 a r 


3.26 a ' 


3.41 a ' 


(|jmol CO2 m' 2 
s- 1 ) 


0.17 


0.25 


0.11 


0.27 


0.13 


0.35 



For each parameter, the mean values ± SE (n = 5-8) followed by different letters 
express significant differences between cultivars for the same CO2 treatment (a, b) 
or between CO2 treatments within the same cultivar (r, s). The ANOVAs for V cmax 



and J max showed significant differences between cultivars for the same CO2 
treatment and between CO2 treatments within the same cultivar. Only for the first 
variable was also detected a significant interaction between genotype and growth 
C0 2 . 

doi: 10.1371/journal.pone.0082712.t002 

Thylakoid electron transport rates 

The potential thylakoid electron transport rate for both 
photosystems changed with growth [C0 2 ], although with 
different extent among the genotypes (Figure 3). Increases in 
the electron transport of photosystem II including the oxygen 
evolving complex (PSII+OEC) under high [C0 2 ] ranged from 
10% (C. arabica genotypes, although this increase was not 
significant in IPR 108) to 30% (CL 153). The PSII activity 
excluding the OEC (PSII-OEC) also showed significant 
increases, ranging from 12% (IPR 108) to 25% (CL 153). The 
pattern of PSI activity closely followed that of PSII, with 
significant enhancements of 20% in CL 153 and 9% in Icatu, 
whereas IPR 108 showed a non-significant rise of 6% in PSI 
activity. 

Enzyme activities 

The enzyme activities were clearly affected by growth [C0 2 ] 
in all genotypes. The potential activities of the photosynthetic 
related enzymes ribulose-1,5-bisphosphate carboxylase/ 
oxygenase (RuBisCo) and ribulose 5-phosphate kinase 



(Ru5PK) were similar among genotypes at normal [COJ. 
Additionally, significant increases were observed for high 
growth [C0 2 ] plants, from 37% (IPR 108) to 46% (Icatu) for 
RuBisCo and from 35% (CL 153) to 63% (IPR 108) for Ru5PK 
(Figure 4). Similar patterns were exhibited by the key enzymes 
of the respiratory pathway, malate dehydrogenase (MDH) and 
pyruvate kinase (PK) (Figure 5), with increases in the plants 
grown under high [COJ ranging from 20% (CL 153) to 75% 
(Icatu) for MDH and from 76% (CL 153) to 86% (Icatu) for PK. 

Non-structural carbohydrate quantification 

Several soluble sugars showed different variation patterns 
under high growth [COJ (Table 4). The most abundant sugar, 
sucrose, increased as much as 20% in IPR 108, whereas 
levels were more stable in CL 153 and Icatu. Glucose did not 
change significantly in CL 153 and Icatu but decreased by 42% 
in IPR 108, whereas fructose was stable in Icatu and 
decreased in the other two genotypes. Raffinose did not 
change significantly in any genotype, and arabinose decreased 
only in CL 153 (32%), in response to higher growth [COJ. 
Finally, strong reductions under high growth [COJ were found 
for trehalose (between 78 and 92%), stachyose (63-88%), 
galactose (58-89%) and mannitol (51-77%). Taken together, 
these changes led to decreases of total soluble sugar content 
in all genotypes, between 15% (Icatu) and 33% (CL 153). 

Starch changes under high [COJ were species-dependent, 
increasing by 69% in CL 153 while decreasing by 10 and 15% 
in Icatu and IPR 108, respectively. 

Total non-structural carbohydrates (NSC) did not increase in 
any of the genotypes. In fact, NSC tended to decrease, 
particularly in IPR 108 (25%) and Icatu (13%). 

Photosynthetic pigments 

Total chlorophylls, total carotenoids and their ratios did not 
differ significantly with respect to genotype or growth C0 2 
treatment (Figure 6). 

Cellular membrane permeability and quantification of 
chloroplast membrane lipids 

Among genotypes, some differences were noted in the 
cellular membrane leakage values at normal growth C0 2 , but 
these were not affected by enhanced growth [COJ (Figure 7). 
In contrast, the lipid fraction displayed significant reductions of 
total fatty acid (TFA) content in CL 153 and Icatu (Figure 8). 
This reduction was accompanied by a decrease in the 
unsaturation level (lower values of the double bond index, DBI) 
in CL 153, although the opposite tendency was found for both 
Icatu and IPR 108. These changes in the DBI of CL 153 
resulted from modifications in the weight of individual FAs 
(Table 5), as the contents of the two most important FAs, 
palmitic acid (C16:0) and linolenic acid (C18:3), followed 
opposite trends, with the latter decreasing significantly under 
high [COJ. In the C. arabica genotypes, the weight of the 
individual FAs remained mostly unaltered, although minor 
changes occurred in C16:0 and C18:3, resulting in small 
increases in DBI. Additionally, the 3-trans-hexadecenoic acid 
(C16:1 c+t) tended to increase, albeit only significantly in IPR 
108. 
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Figure 2. Potential photosynthesis. Changes in photosynthetic capacity (A max ) in the leaves of C. arabica (Icatu and IPR 108) 
and C. canephora (Conilon CL 153) grown under 380 (white bar) and 700 (black bar) uL C0 2 L 1 . The mean values + SE (n = 6) 
followed by different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 
treatments within the same cultivar (r, s). The ANOVA for A max showed significant differences between cultivars for the same growth 
C0 2 treatment. 

doi: 10.1371/journal.pone.0082712.g002 



Discussion 

Long-term [C0 2 ] enrichment did not provoke down- 
regulation of either stomatal conductance or 
photosynthesis 

Stomatal size and density are key determinants of maximum 
g s [38,54]. Here, we showed that the dichotomous behavior of 
these traits likely contributed to the low response of g s to [C0 2 ] 
enrichment in the tested genotypes (Table 1). These results 
agree with those of few studies that have reported invariant 
responses of g s to elevated [C0 2 ] in woody species [8,55,56]. 
Indeed, g s has been systematically, but not universally, 
demonstrated to decrease when plants are grown under 
elevated [C0 2 ] [8,16,17,22], although such decreases have 
been shown to occur to a lesser extent in shrubs and trees 
than in herbaceous annuals [8]. 

In the current study, no apparent photosynthetic down- 
regulation to long-term elevated [C0 2 ] was observed. Rather, 
some components of the photosynthetic machinery were even 
up-regulated under high [COJ. The present work presents 
compelling evidence that supports this conclusion. First, 
significant increases in P n were found in either genotype when 
measured under elevated [C0 2 ] (Figure 1). These increases 
were likely related not only to a higher carboxylation rate linked 



Table 3. Evaluation of several fluorescence parameters 
determined under dark adapted (F V /FJ and under 
photosynthetic steady-state (q P , F v 7F m ') conditions in the 
leaves of C. arabica (Icatu and IPR 108) and C. canephora 
(Conilon CL 153) grown under 380 and 700 uLC0 2 L 1 . 



Genotype 


CL 153 




Icatu 




IPR 108 




Growth 


380 ul_ 


700 uL 


380 uL 


700 uL 


380 uL 


700 uL 


[C0 2 ] 


L' 1 


L- 1 


L- 1 


L" 1 


L- 1 


L- 1 


F v /F m 


0.776 a ' 


0.760 a ' 


0.754 b ' 


0.732 b s 


0.758 b ' 


0.748 ab ' 




0.006 


0.007 


0.006 


0.006 


0.004 


0.005 


qp 


0.682 ar 


0.558 a ' 


0.659 a ' 


0.572 a ' 


0.590 a ' 


0.621 a ' 




0.022 


0.026 


0.060 


0.069 


0.024 


0.034 


F v ''F m ' 


0.575 ar 


0.577 ab ' 


0.589 a ' 


0.664 a ' 


0.531 ar 


0.542 b r 




0.015 


0.024 


0.028 


0.024 


0.026 


0.023 



For each parameter, the mean values ± SE (n = 5) followed by different letters 
express significant differences between cultivars for the same CO2 treatment (a, b) 
or between C0 2 treatments within the same cultivar (r, s). Only the ANOVA for qp 
showed significant effects between C0 2 treatments within the same cultivar. 
doi: 10.1371/journal. pone. 008271 2.t003 
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Figure 3. Thylakoid electron transport potential. Changes in the thylakoid electron transport rates associated with PSII (with 
and without the inclusion of OEC) and PSI in the leaves of C. arabica (Icatu and IPR 108) and C. canephora (Conilon CL 153) 
grown under 380 (white bar) and 700 (black bar) uL C0 2 L 1 . For each parameter, the mean values + SE (n = 4) followed by different 
letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 treatments within the 
same cultivar (r, s). The ANOVAs for PSII+OEC, PSII-OEC and PSI showed significant differences between cultivars for the same 
C0 2 treatment and between C0 2 treatments within the same cultivar, as well for the interaction between genotype and growth C0 2 
(except in PSI for the interaction). 

doi: 10.1371/joumal.pone.0082712.g003 
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Figure 4. Maximal activities of some photosynthetic enzymes. Variation of the total activities of ribulose-1 ,5-bisphosphate 
carboxylase/oxygenase (RuBisCo) and ribulose 5-phosphate kinase (Ru5PK) in the leaves of C. arabica (Icatu and IPR 108) and C. 
canephora (Conilon CL 153) grown under 380 (white bar) and 700 (black bar) uL C0 2 L" 1 . For each enzyme, the mean values + SE 
(n = 4) followed by different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between 
C0 2 treatments within the same cultivar (r, s). The ANOVA for RuBisCo showed significant differences only between C0 2 
treatments within the same cultivar, whereas that for Ru5PK showed significant differences between cultivars for the same C0 2 
treatment and between C0 2 treatments within the same cultivar. 

doi: 10.1371/joumal.pone.0082712.g004 
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Figure 5. Maximal activities of some respiratory enzymes. Variation of the total activities of malate dehydrogenase (MDH) and 
pyruvate kinase (PK) in the leaves of C. arabica (Icatu and IPR 108) and C. canephora (Conilon CL 153) grown under 380 (white 
bar) and 700 (black bar) uL C0 2 L 1 . For each enzyme, the mean values + SE (n = 4) followed by different letters express significant 
differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 treatments within the same cultivar (r, s). The 
ANOVAs for MDH and PK showed significant differences between the cultivars for the same C0 2 treatment and between C0 2 
treatments within the same cultivar. 

doi: 10.1371/joumal.pone.0082712.g005 
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Table 4. The concentrations of non-structural 
carbohydrates (soluble sugars and starch) in the leaves of 
C. arabica (Icatu and IPR 108) and C. canephora (Conilon 
CL 153) grown under 380 or 700 uL C0 2 L 1 . 



Genotype 


CL 153 




Icatu 




IPR 108 




co 2 


380 uL 


700 uL 


380 uL 


700 uL 


380 uL 


700 uL 


Treatment 


L" 1 


L- 1 


L" 1 


L" 1 


L- 1 


L" 1 


Stachyose 
(mg g" 1 dw) 


1.35 b r 


0.28 ab s 


1.11 c ' 


0.41 a s 


1.60 a r 


0.20 b r 




0.04 


0.02 


0.08 


0.01 


0.03 


0.01 


Raffinose (mg 

g" 1 dw) 


3.54 c r 


3.65 b r 


4.06 b r 


3.62 b r 


4.98 a r 


5.11 a r 




0.12 


0.12 


0.13 


0.18 


0.08 


0.10 


Trehalose 
(mg g -1 dw) 


8.12 ab r 


1 .09 ab s 


7.61 b r 


1.69 a 5 


8.64 a r 


0.73 b 5 




0.24 


0.08 


0.54 


0.03 


0.15 


0.04 


Sucrose (mg 

g" 1 dw) 


38.76 b r 


37.67 b r 


44.12 a ' 


41.52 a r 


32.12 c s 


38.70 ab 
r 




1.06 


0.35 


0.78 


1.68 


0.93 


0.39 


Glucose (mg 

g" 1 dw) 


4.69 c ' 


5.23 c r 


16.87 b r 


16.92 a r 


25.53 a r 


14.80 a 5 




0.24 


0.19 


0.40 


0.44 


1.01 


0.40 


Fructose (mg 

g" 1 dw) 


23.31 bc 
r 


11.47 b 5 


21.57 c r 


19.91 a r 


26.93 a r 


17.60 a s 




0.45 


0.38 


0.59 


1.45 


1.53 


0.40 


Galactose 
(mg g -1 dw) 


0.47 b r 


0.05 b 5 


0.38 c ' 


0.16 as 


0.72 a ' 


0.12 a 5 




0.01 


0.00 


0.03 


0.00 


0.01 


0.01 


Arabinose 
(mg g" 1 dw) 


4.74 a ' 


3.19 a s 


4.38 ab r 


3.88 a r 


3.76 b r 


3.13 a r 




0.07 


0.28 


0.06 


0.34 


0.31 


0.24 


Mannitol (mg 
g" 1 dw) 


11.23 a r 


2.58 b s 


9.67 b r 


4.73 a s 


10.70 ab 
r 


2.59 b s 




0.33 


0.20 


0.69 


0.08 


0.19 


0.16 


Total Soluble 
(mg g" 1 dw) 


96.20 b ' 


64.86 c s 


109.77 a 
r 


92.83 a 5 


114.99 a 
r 


82.98 b 5 




1.60 


0.54 


2.57 


0.75 


2.83 


0.36 


Starch (mg 














glue. 

equivalents 9 ^ 


33.16 b 5 


56.32 a ' 


54.68 a ' 


49.21 a ' 


25.09 b r 


21.44 br 














dw) 
















1.80 


5.73 


1.53 


5.07 


0.81 


4.20 


Total Soluble/ 














Starch (g g" 1 ) 


2.90 


1.15 


2.01 


1.89 


4.58 


3.87 


Total Sugars 
(mg g -1 dw) * 


129.36 


121.18 


164.45 


142.04 


140.08 


104.42 



to the increase of C0 2 as a substrate but also to the 
competitive inhibition of the oxygenation reaction of RuBisCo 
and the subsequent reduction of C0 2 loss and energy costs 
associated with the photo respiratory pathway [7,8]. Second, 
despite the increase in P n , at least a partial photosynthetic and 
stomatal down-regulation could have occurred due to 
prolonged growth under high [COJ. In fact, this was not the 
case; the P n increase was independent of growth [C0 2 ], as no 



Table 4 (continued). 



For each sugar, the mean values ± SE (n = 6) followed by different letters express 
significant differences between cultivars for the same C0 2 treatment (a, b, c) or 
between C0 2 treatments within the same cultivar (r, s). The ANOVAs for 
stachyose, glucose, fructose and galactose showed significant differences 
between cultivars for the same C0 2 treatment, between C0 2 treatments within the 
same cultivar, and for the interaction between genotype and growth C0 2 ; that for 
arabinose showed significant differences between cultivars for the same C0 2 
treatment and between C0 2 treatments within the same cultivar; those for sucrose 
and starch showed significant differences between cultivars for the same C0 2 
treatment and for the interaction; those for trehalose and mannitol showed 
significant differences between C0 2 treatments within the same cultivar and for the 
interaction; and that for raffinose showed significant differences only between 
cultivars for the same C0 2 treatment. 

* Values obtained using the mean values of each sugar compound, 
doi: 10.1 371/journal. pone. 0082712.t004 

significant differences existed between plants grown under 
normal or high [C0 2 ] when measurements were performed at 
380 or 700 pL C0 2 L 1 (data not shown). This result indicates 
that the P n stimulation was similar to the potential maximum 
stimulation obtained by the modeled P n /Ci response under 
current [COJ. Third, A max (assessed under saturating C0 2 in 
the absence of diffusion-mediated limitations of 
photosynthesis) was unresponsive to the [C0 2 ] treatments 
(Figure 2). Fourth, the acclimation of the key photosynthetic 
parameters (V cmax and J ma)( ) was not observed (Table 2), in 
sharp contrast to previous studies [for review, see 4,14]. In fact, 
photosynthesis is known to shift from limitation by Rubisco 
carboxylation capacity to consume RuBP at low Ci (because 
C0 2 is the substrate), to RuBP regeneration-limited rates at 
higher [C0 2 ], related to the ability of thylakoid-dependent 
reactions to supply ATP and NAPDH [8,10], as well as to the 
capacity of starch and sucrose synthesis to utilize those 
phosphates and subsequently regenerate P, for 
photophosphorylation [21]. Irrespective of these assumptions, a 
modest average decrease (6%) of V cmax has been observed in 
trees under high [C0 2 ] [8], whereas a downward adjustment of 
J max has been found less frequently [7]. The absence of the 
acclimation of J max and V cmax [9], or a positive acclimation to 
elevated [C0 2 ], could also occur, shifting the investment of 
resources from RuBisCo to the processes supporting RuBP 
and Pi regeneration. Fifth, TPU, which may become limited 
under high [COJ even without photosynthetic acclimation in 
trees [9], remained unchanged. This result suggests that P n is 
not limited by the plant's ability to synthesize (and use) starch 
and sucrose, using those phosphates and regenerate P; [21]. 
Sixth, the thylakoid electron transport involving both PSII and 
PSI, as well as the photosynthetic enzyme activities, were, 
overall, up-regulated under elevated [COJ (Figure 3). Finally, 
no apparent up-regulation of the NSC pools was observed (see 
below). Considered together, these results suggest that if root 
growth is unrestricted by pot size [26,57], as it was in the 
present case, the coffee plant can sustain relatively high 
photosynthetic rates when growing under elevated [C0 2 ] 
whenever water, temperature and nutrients are not limiting. 
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Figure 6. Photosynthetic pigment stability. Changes in the content of total chlorophylls (Chi a+b) and total carotenoids (Tot. 
Carot.), as well as their ratio (Tot Chi/Tot Carot), in the leaves of C. arabica (Icatu and IPR 108) and C. canephora (Conilon CL 153) 
grown under 380 (white bar) and 700 (black bar) uL C0 2 L 1 . For each parameter, the mean values + SE (n = 6-8) followed by 
different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 treatments 
within the same cultivar (r, s). 

doi: 10.1371/journal.pone.0082712.g006 
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Figure 7. Cellular membranes selectivity stability. Evaluation of membrane permeability in the leaves of C. arabica (Icatu and 
IPR 108) and C. canephora (Conilon CL 153) grown under 380 (white bar) and 700 (black bar) uL C0 2 L-\ The mean values + SE (n 
= 5) followed by different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 
treatments within the same cultivar (r, s). The ANOVA for leakage showed significant differences between cultivars for the same 
C0 2 treatment, between C0 2 treatments within the same cultivar and for the interaction. 

doi: 10.1371/joumal.pone.0082712.g007 



Another widely reported effect promoted by high growth 
[C0 2 ] is the increase of WUE, often resulting from the 
maintenance of or decreases in g s (and T r ) in parallel to 
increases in P n , even when partial photosynthetic down- 
regulation occurs [8,14,17,19]. In this study, the increase in 
iWUE in CL 153 was mostly promoted by the enhancement of 
P n , as has been observed in grapevine [56] and in forest trees 
[9], whereas in Icatu and IPR 108 a tendency to lower g s values 
also contributed to iWUE increase (Figure 1 ). Since changes in 
precipitation patterns and an increased frequency of drought 
episodes are predicted due to global warming [1,7], the 
improved leaf-level iWUE of plants grown under enhanced 
[C0 2 ] could be of great interest, as it may help protect plants 
from water shortage. However, it must be noted that WUE 
increases due to transpiration (and g s ) reductions at the leaf 
level do not necessarily lead to corresponding reductions in 
conductance and evapotranspiration at the canopy level [58]. 
Moreover, in the CL 153 plants, iWUE increased mostly due to 
a P n upsurge under high C0 2 conditions. Knowing that the 
predicted rise in temperature (and lower rainfall) that will 
accompany the increase in atmospheric C0 2 , is expected to 
increase the leaf-to-air vapor pressure (that is a crucial 
determinant of leaf transpiration, together with g s ), thus, no 
large reduction of water use would be expected in CL 153, 



contrary to what could happens in the C. arabica genotypes 
where g s tend to lower values. 

Metabolic machinery was adjusted in response to high 
[COJ 

In all tested genotypes, the most prominent effect of [C0 2 ] 
was the up-regulation of the in vitro activities of key enzymes of 
carbon metabolism, Ru5PK (a crucial enzyme in the RuBP 
regeneration pathway), RuBisCo (Figure 4), MDH and PK 
(Figure 5), suggesting a reinforcement of the potential 
biochemical capacities of both photosynthesis and respiration 
under the present experimental conditions. Such reinforcement 
might prevent the down-regulation of photosynthesis, which 
has been associated with reductions in N allocation to 
RuBisCo, RuBP regeneration and proteins associated with 
electron transport [9]. However, the increase of total RuBisCo 
activity did not match the behavior of V cmax in CL 153, which 
might be explained by a decreased RuBisCo activation state. 
Accordingly, increased activities of MDH and PK were also not 
accompanied by increases in respiration rates (R d ), which 
remained unchanged in response to high [COJ (Figure 2). 
Indeed, depending on the species, declining or unchanging R d 
have commonly been observed under high [COJ [see 3,16,59]. 

Increases in P n under high [COJ obviously require more 
energy and reducing equivalents, which is consistent with the 
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Figure 8. Dynamics of the lipid matrix from chloroplast membranes. Changes in the total fatty acid (TFA) content and 
unsaturation level (DBI) of the chloroplast membranes of leaves of C. arabica (Icatu and IPR 108) and C. canephora (Conilon CL 
153) grown under 380 (white bar) and 700 (black bar) uL C0 2 L 1 . For each fatty acid, the mean values + SE (n = 4) followed by 
different letters express significant differences between cultivars for the same C0 2 treatment (a, b) or between C0 2 treatments 
within the same cultivar (r, s). The ANOVAs for TFA and DBI showed significant differences between cultivars for the same C0 2 
treatment, between C0 2 treatments within the same cultivar, and for the interaction between genotype and growth C0 2 . 

doi: 10.1371/joumal.pone.0082712.g008 
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Table 5. Evaluation of the proportions of palmitic acid 
(C16:0), 3-trans-hexadecenoic acid (C16:1t), stearic acid 
(C18:0), oleic acid (C18:1), linoleic acid (C18:2) and 
linolenic acid (C18:3), from the chloroplast membranes of 
the leaves of C. arabica (Icatu and IPR 108) and C. 
canephora (Conilon CL 153) grown under 380 and 700 uL 
C0 2 L- 1 . 



Genotype 


CL 153 




Icatu 




IPR 108 




co 2 


380 |JL 


700 uL 


380 uL 


700 uL 


380 uL 


700 uL 


Treatment 


L" 1 


L' 1 


L" 1 


L- 1 


L- 1 


L" 1 


<C16:0 (%) 


2.54 c s 


5.62 b ' 


4.57 b 5 


7.56 a r 


5.33 a ' 


4.80 c ' 




0.32 


0.10 


0.32 


0.08 


0.31 


0.18 


C16:0 (%) 


21.9 cr 


24.6 b ' 


32.3 a ' 


27.8 a s 


26.8 br 


24.8 b ' 




1.2 


0.2 


1.1 


1.1 


1.3 


0.6 


C16:1c+t (%) 


3.26 a ' 


2.10 br 


2.31 ab ' 


4.07 a ' 


1.69 bs 


4.43 a ' 




0.55 


0.81 


0.28 


0.34 


0.38 


0.51 


C18:0 (%) 


7.42 a r 


5.86 a r 


7.88 a ' 


7.18 ar 


7.53 a ' 


6.80 a ' 




0.27 


0.19 


0.99 


0.52 


0.11 


0.29 


C18:1c+t (%) 


1.71 bs 


2.13 ar 


2.30 a ' 


1.95 ar 


2.26 ar 


1.59 b 5 




0.01 


0.02 


0.18 


0.16 


0.03 


0.03 


C18:2 (%) 


11.1 bs 


14.1 ar 


14.2 a ' 


12.8 br 


13.7 ar 


13.5 abr 




0.2 


0.3 


0.8 


0.8 


0.4 


0.4 


C18:3 (%) 


52.1 a ' 


45.6 a s 


36.4 c ' 


38.6 br 


42.7 b ' 


44.1 ar 




1.2 


0.9 


2.0 


0.1 


0.9 


0.3 



For each fatty acid, the mean values ± SE (n = 4) followed by different letters 
express significant differences between cultivars for the same CO2 treatment (a, b, 
c) or between CO2 treatments within the same cultivar (r, s). The ANOVAs for all 
fatty acids showed significant differences between cultivars for the same CO2 
treatment, between CO2 treatments within the same cultivar and for the interaction 
between genotype and growth CO2, except that for C18:0 that did not present a 
significant interaction effect, 
doi: 10.1371/joumal.pone.0082712.t005 

increases in the potential electron transport capacity (Figure 3). 
These parallel increases were likely associated with the 
maintenance of a close functional balance between 
carboxylation and electron transport events (J max /V cmax ), which 
seemed to be conserved among coffee genotypes and 
irrespective of growth [C0 2 ] (Table 2). Similar results have 
been reported for some acacia species [17,60] and also in 
other plant species, where the unchanging J max /V cmax ratio has 
been interpreted as reflecting an absence of resource 
redistribution among photosynthetic components [61]. 
Additionally, the reinforcement of photosynthetic components 
under high [C0 2 ] was also unrelated to changes in 
photochemical efficiency as analyzed by Chi a fluorescence. Of 
particular interest, F v /F m , F v 7F m ' and q P (Table 3) were quite 
stable across experimental conditions, as has also been found 
in grapevine [56]. Notably, these results agree with the 
absence of noticeable changes in the pools of total chlorophyll 
(Chi) and carotenoids (Figure 6), as has also been reported for 
Chi in other woody species [7,9]. Altogether, these data 
indicate a pattern of enhanced investment in key components 
of both photosynthetic and respiratory pathways under high 
growth [C0 2 ] in coffee plants. 



Absence of non-structural sugars accumulation 

The accumulation of non-structural carbohydrates (NSC) in 
leaves is one of the most pronounced and common responses 
of C 3 plants to elevated [C0 2 ], even in field-grown plants where 
rooting volume is unrestricted [7]. In coffee genotypes, several 
soluble sugars presented significant reductions under high 
growth [C0 2 ] conditions, particularly trehalose, stachyose, 
galactose and mannitol, whereas fructose and glucose either 
remained stable or decreased, resulting in a reduction of both 
total soluble sugars and NSC contents (Table 4). Therefore, 
the absence of sugar accumulation in leaves would have 
contributed to the avoidance of photosynthetic down-regulation 
(Figure 1, Table 2), and further supports the lack of negative 
impacts on the photosynthetic enzymes studied (Figure 4). 
Similar NSC findings have been reported in young sunflower 
leaves [24], as well as in poplar, which was able to export over 
90% of its photosynthate during the day and had a large 
capacity for temporary storage in starch, thus, maintaining its 
potential for C acquisition [62]. 

In coffee, the absence of an NSC increase could be related 
to the higher production of new leaves and plagiotropic (lateral) 
branches (somewhat modifying the architecture of the coffee 
plant), as well as flower and fruit production occurring up to 
twice a year (unpublished observation). This higher continuous 
production of vegetative and reproductive structures would 
maintain a higher consumption of photosynthates, contributing 
to an increased sink strength that would ultimately explain the 
absence of photosynthetic acclimation. Enhanced sink strength 
linked to the maintenance of the photosynthetic stimulation 
under C0 2 enrichment has been observed in other species 
such as Vernonia herbacea [6], loblolly pine [63], sour orange 
[4] and grapevine [56], with positive yield implications. 

Preservation of cell membrane permeability and 
changes in chloroplast lipids 

Under normal [C0 2 ], TFA content and the degree of 
unsaturation (DBI) were slightly higher in the C. canephora 
genotype (Figure 8), as has also been previously reported [28]. 
Cell membrane permeability was not significantly modified due 
to growth [C0 2 ] (Figure 7), but the lipid matrix of the chloroplast 
membranes was altered in terms of TFA content in CL 153 and 
Icatu. Additionally, variations in the major C16:0 and C18:3 
FAs, which increased and decreased, respectively (Table 5), 
led to a decline in the lipid membrane unsaturation in CL 153, 
whereas the opposite tendency was found in the C. arabica 
genotypes, with possible implications for membrane fluidity 
[28,64]. Furthermore, C1 6: 1c+f values tended to increase in the 
C. arabica genotypes, becoming somewhat higher than in CL 
153. The C16:1t is a major and specific FA of the 
phosphatidylglycerol (PG) in chloroplasts [28,65,66]. Both PG 
and C16:1t contribute to preserve the thylakoid membrane's 
organization of proteins and pigments and the optimal 
conformation of the D1 protein, stabilizing the photosystem 
complexes and allowing for efficient non-cyclic electron flow 
[66,67]. However, although some differences among genotypes 
were observed under high [C0 2 ] concerning the lipid 
unsaturation degree and the composition in some FAs in the 
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chloroplast membranes, no distinct impacts on plant function, 
e.g., changes in electron transport, could be detected. 

No clear differences between C. arabica and C. 
Canephora genotypes were observed in response to 
elevated [C0 2 ] 

No clear species-dependent responses to elevated [C0 2 ] 
were found, as the three genotypes displayed similar trends for 
most of the collected data. Nonetheless, some differences 
seemed to exist between the genotypes. For example, IPR 108 
was somewhat less responsive for most parameters, Icatu 
presented the highest increases on V cmax and J max , whereas CL 
153 showed the more stable V cmax and J max values (Table 2), 
the highest thylakoid electron transport capability (Figure 3) 
and starch increase (Table 4), accompanied by a reduction in 
the unsaturation of chloroplast membrane lipids (Figure 8, 
Table 5). Regardless of these differences and genotypes, our 
data suggest that coffee can successfully cope with high [C0 2 ] 
under optimal conditions (water, temperature and nutrient 
availability), but further studies are necessary to evaluate its 
responses to a broader scenario of climatic changes, including 
altered water availability and high temperatures. 

Conclusions 

Under enhanced growth [C0 2 ], no down regulation of g s and 
P n was observed in the current study under unrestricted 
conditions of water, nutrients and root development. Stomatal 
density and size showed dichotomous behaviors, decreasing 
and increasing, respectively, whereas g s exhibit no significant 
response to elevated [COJ. Most importantly, changes in P n 
largely governed the significant rise in the iWUE under 
elevated [C0 2 ], although in the C. arabica cultivars g s may have 
somewhat contributed. No negative impacts on the P n , J max and 
Vcmax values were observed, a result most likely linked with the 
absence of NSC accumulation. Considering the stable 
environmental conditions of this experiment, this absence of 
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